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Abstract 

A (n + l)-dimensional Einstein-Gauss-Bonnet (EGB) model is considered. For 
diagonal cosmological-type metrics, the equations of motion are reduced to a set of 
Lagrange equations. The effective Lagrangian contains two "minisuperspace" met- 
rics on W 1 . The first one is the well-known 2- metric of pseudo-Euclidean signature 
and the second one is the Finslerian 4-metric that is proportional to n-dimensional 
Berwald-Moor 4-metric. When a "synchronous-like" time gauge is considered the 
equations of motion are reduced to an autonomous system of first-order differen- 
tial equations. For the case of the "pure" Gauss-Bonnet model, two exact solu- 
tions with power-law and exponential dependence of scale factors (with respect to 
"synchronous-like" variable) are obtained. (In the cosmological case the power-law 
solution was considered earlier in papers of N. Deruelle, A. Toporensky, P. Tretyakov 
and S. Pavluchenko.) A generalization of the effective Lagrangian to the Lowelock 
case is conjectured. This hypothesis implies existence of exact solutions with power- 
law and exponential dependence of scale factors for the "pure" Lowelock model of 
m-th order. 
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1 Introduction 



Here we deal with D-dimensional gravitational model with the Gauss-Bonnet term. The 
action reads 



S= [ d D z^\g~\{ ai R[g] + a 2 C 2 [g}}, (1.1) 

J M 



where g = guNdz M ® dz N is the metric defined on the manifold M, dimM = D, \g\ = 
I det(g MN ) \ and 

C 2 = RmnpqR MNPQ — 4:RmnR MN + R 2 (1-2) 

is the standard Gauss-Bonnet term. Here a± and a 2 are constants. The appearance of 
the renormalizable Gauss-Bonnet term as well as quadratic Riemann curvature terms in 
multidimensional gravity is motivated by string theory jTJ |2j [3j HJ 15] ■ (For a review of 
fourth-order gravity in D = 4, see [6].) 

At present, the so-called Einstein-Gauss-Bonnet (EGB) gravity and its modifications 
are intensively used in cosmology, see [S] (for D = 4), [EH CE21 EH (H H51 US] 
and references therein, e.g. for explanation of accelerating expansion of the Universe 
following from supernovae (type la) observational data [17] . Certain exact solutions in 
multidimesional EGB cosmology were obtained in [9]- [16] and some other papers. 

EGB gravity is also intensively investigated in a context of black-hole physics. The 
most important results here are related with the well-known Boulware-Deser- Wheeler 
solution (corresponding to the Schwarzschild-Tangherlini solution in general relativity) 
[T8| [T9] and its generalizations [201 IZH [22j [23] , for a review and references, see [2U [25] . 
For certain applications of brane-world models with Gauss-Bonnet term, see review [26] 
and references therein. 

Here we are interested in the cosmological (type) solutions with diagonal metrics (of 
Bianchi-I-like type) governed by scale factors depending upon one variable. 

For a 2 = 0, we have the Kasner type solution with the metric 

n 

g = -dT(g)dT + J2 4 2 t V cV ® dy\ (1.3) 

8=1 

where Ai > are arbitrary constants, D = n + 1 and parameters p l obey the relations 

n 

$> = 1, (1.4) 

2=1 

n 
i=l 

and hence 

n i n 

E P i P i= 2^P^-2^ {pi)2 = °- (L6) 

l<i<j'<n i=l i=l 
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For D = 4, this is the well-known Kasner solution [27]. The set of eqs. ( II. 4p . ( II. 5p is 
equivalent to the set of eqs. ( II. 4p . ( II. 6p . 

In [10], a Einstein- Gauss-Bonnet (EGB) cosmological model was considered. For 
"pure" Gauss-Bonnet (GB) case a.\ = and c*2 ^ 0, N. Deruelle has obtained a cos- 
mological solution with the metric (II. 3p for n = 4, 5 and parameters obeying the relations 



£y=3, (i.7) 



i=l 



pYp k p l = 0. (1.8) 

l<i<j<k<l<n 

It was reported by A. Toporensky and P. Tretyakov in [13j that this solution was verified 
by them for n = 6, 7. In recent paper by S. Pavluchenko [28J the power-law solution was 
verified for all n (and also generalized to the Lowelock case [3"4]). 

In this paper we give a derivation of the "power-law" (cosmological type) solution for 
arbitrary n. We also show that for D ^ 4 this solution in "pure" GB cosmology is unique 
in a class of solutions with power-law dependence of scale factors: aj(r) = AiT p \ when the 
parameters p 1 , ...,p n contain more than two non-zero numbers. When (n — 2) parameters 
among p % are zero, say p 3 = ... = p n = 0, than the metric (II. 3p obeys the equations of 
motion (for ot\ = 0) for arbitrary values of two Kasner-like parameters (say p l ,p 2 ). 

The numerical analysis of cosmological solutions in EGB gravity for D = 5, 6 [15] shows 
that the singular "power-law" solutions (II. 3p . ( II. 7p . (11. 8p (e.g. with a little generalization 
of scale factors a,i{r) = Aj(r ±r) p \ where r is constant) appear as asymptotical solutions 
for certain initial values as well as Kasner- type solutions ( I1.3p -( 1T31) do. 

The paper is organized as follows. In Section 2 the equations of motion for (n + 1)- 
dimensional EGB model are considered. For diagonal cosmological type metrics the equa- 
tions of motion are reduced to a set of Lagrange equations corresponding to certain "effec- 
tive" Lagrangian (in agreement with [iTH [28] for cosmological case). Section 3 is devoted 
to the case of the "pure" Gauss-Bonnet model. Two exact solutions: with power-law and 
exponential dependence of scale factors (with respect to "synchronous-like" variable) are 
obtained. In Section 4 the equations of motion are reduced to an autonomous system 
of first order differential equations (when a "synchronous-like" time gauge is considered). 
For «! 7^ and a<i ^ it is shown that for any non-trivial solution with exponential 
dependence of scale factors cl^t) = Ajexp(fV), i = l,...,n, there are no more than 
three different numbers among v 1 ,...,v n . In Section 5 a generalization of the effective 
Lagrangian to the Lowelock case is conjectured and exact solutions with power-law and 
exponential dependence of scale factors for the "pure" Lowelock model of m-th order 
are presented. (See also [10], [28] for "power law" cosmological solutions.) Certain useful 
relations and proofs are collected in Appendix. 
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2 The cosmological type model and its effective La- 
grangian 

2.1 The set-up 

Here we consider the manifold 

M = x Mi x . . . x M n , (2.1) 

with the metric 

n 

g = we 2 ^du ®du + J2 e 2 ^ {u) eidy l ® dy\ (2.2) 

t=i 

where w — ±1 and any Mj is 1-dimensional manifold with the metric g % = £idy l ® dy l , 
Si = ±1, % = 1, . . . ,n. Here and in what follows IR* = (u-,u + ) is an open subset in R. 
(Here we identify g' with g l = p*g l which is the pullback of the metric g l to the manifold 
M by the canonical projection: : M — > Mj, i = l,...,n.) The functions j(u) and 
/3 l (u), i = 1, . . . , n, are smooth on = (w_, 

For w = —1, Ex = ... = e n = 1 the metric ( 12. 2 p is a cosmological one while for w = 1, 
£i = — 1, e 2 = ... = E n = 1 it describes static configurations. 

According to Appendix A, the integrand in (II. ip . when the metric ( 12. 2 p is substituted, 
reads as follows 

y/\a\{ ai R\g] + a 2 C 2 [g]} = L + |-, (2.3) 



where 



du 



L = aiLi + a 2 L 2 , (2.4) 
Li = (-u0e-^a tJ /3^, (2.5) 

= -\e-^G i3kl frfcp k p\ (2.6) 



7o = EILi P and 



Gij = 5ij - 1, (2.7) 
G m = {5 l3 - l)(5 ik - l)(5 a - l)(5 jk - 1)(^ - l)^ - 1) (2.8) 

are respectively the components of two "minisuperspace" metrics on M. n . (For cosmo- 
logical case see also pUl EHl E7].) The first one is the well-known 2-metric of pseudo- 
Euclidean signature: < vi,v 2 >= Gijv\v\ and the second one is the Finslerian 4-metric: 
< t>i, t>2, t>3, t>4 >= Gij k iv\v 2 v\v\, v s = (v l s ) G M n , where < .,. > and < > are 

respectively 2- and 4-linear symmetric forms on IR n . (Here we denote A = dA/du etc.) 
In (12. 3p the function / = /(7, (3, $) has the following form: 

f = aifi + a 2 f 2 , (2.9) 
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where fi and / 2 are defined in Appendix A (see (IA.22|) and (1A.23|) ). 

The derivation of (j2.4p - (l'2.6p is based on the relations obtained in Appendix A (see 
f EOU]) . ipijl ) and the following identities 



n 



G^V=^(^) 2 -E^) 2 , (2-10) 



i=l i=l 
n n 



j=i i=i j=i 

n n n n 

+3£>?) 2 + 8(JV) J> J ) 3 - 6^K) 4 - (2.11) 

i=l i=l j=l i=l 

The first identity (I2.10p is a trivial one. The second one (12. lip may be verified by straight- 
forward calculations (see Appendix B). 

It follows immediately from the definitions ( 12. 7p and (12. 8p that 

GijvW = -2^uV, (2.12) 

i<j 

G ljk iv i v j v k v l =24 ^ «Vt)V. (2.13) 

i<j<fc<Z 

Due to (I2~T31 . G^bVW is zero for n = 1,2,3 (D = 2,3,4). For n = A (D = 5), 
Gijuv % v^v k v l = 24v l v 2 v 3 v 4 and our 4-metric is proportional to the well-known Berwald- 
Moor 4-metric [291 EQ] (see also [311 [32] and references therein). We remind the reader 
that the 4-dimensional Berwald-Moor 4-metric obeys the relation: < v,v,v,v >bm= 
v 1 v 2 v 3 v A . The Finslerian 4-metric with components ( 12. 8 p coincides up to a factor with 
the n-dimensional analogue of the Berwald-Moor 4-metric. 



2.2 The equations of motion 

The equations of motion corresponding to the action (II. ip have the following form 

S MN = a x S$ N + a 2 S$ N = 0, (2.14) 

where 

&mn = Rmn — -RgMN, (2.15) 

c( 2 ) _ Of T? E> P Q S Of? E> P 

—2R MPN qR pq + RR MN ) — -C 2 gMN- (2.16) 

The field equations (I2.14p for the metric (12.21) are equivalent to the Lagrange equations 
corresponding to the Lagrangian L from (12.41) . This follows from the relations 

dL 

S 00 (-2w) exp( 7o - 7 ) = — , (2.17) 



£jj (- 2£i )exp( 7 + 70 -2« = ^--_, (2.18, 

£m = 0, (2.19) 

t = 1, . . . ,n. 

Formulas fl2.17p - fl2.19p may be verified just by straightforward calculations based on 
the relations for the Riemann tensor from Appendix A. But there exists a more "economic" 
way to prove these formulas using: (i) the diagonality of the matrix Emn (in coordinates 
(y M ) = (y° = u iV % y)'-, (ii) t ne dependence of this matrix only on one variable it, i.e. 
£mn = £mn(u); (iii) the relation ( 12. 31) . The proof of ( 12.17p -( !2~T9l) is given in Appendix 
C. 

Thus, equations (I2.14p read as follows 

wa 1 G ij $ i $ j + otze-^Gijufrft^p 1 = 0, (2.20) 

-^-[-2w« 1 G i7 -e- 7+70 /3 J 
du 

~a 2 e-^G ljkl frp k p l ] - L = 0, (2.21) 

i = 1, . . . , n. Due to fl2T20l 

2 



-w 



-e-^atGijFF. (2.22) 



3 Exact solutions in Gauss-Bonnet model 

Now we put cci = and a 2 7^ 0, i.e. we consider the cosmological type model governed 
by the action 

S 2 = a 2 [ d D z^\g~\C 2 [g}. (3.1) 

J M 

The equations of motion (12. 14ft in this case read 



j (2) _-o(2) 

where 



&mn — K-mn ~ 2^ 2 9mn — 0, (3.2) 



^MiV ~~ 2(RmpqsRn PQS ~ ^RmpRn P 

—2R M pnqR PQ + RRmn)- (3.3) 

Due to identity g MN 1l$ N = 2£ 2 , the set of eqs. (JOJ for D ^ 4 implies 

£2 = 0. (3.4) 

It is obvious that the set of eqs. (13.21) is equivalent for D 7^ 4 to the following set of 
equations 

K% = 0. (3.5) 
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Equations of motion ( I2.20p and (I2.2ip in this case read as follows 



Gijkifrft P k ft = 0, (3.6) 

d 



du . 



e 



37+70 G t]kl frp k $ l 



0, (3.7) 



i — 1, . . . , n. Here L = due to (13.61) . 

Let us put (3 l = for all % or, equivalently, 



= du + 4, (3.8) 
where d and d are constants, i = 1, . . . , n. We also put 



37 = 7o = X> i > ( 3 - 9 ) 
i=i 

i.e. a modified "harmonic" variable is used. Recall that in the case ot\ ^ and «2 = 0, 
the choice 7 = 70 corresponds to the harmonic variable u [33] . 

Then eqs. (13. 7p are satisfied identically and eq. (13. 6p gives us the following constraint 

G ijk ic i c j c k c l = 24 cVcV = 0. (3.10) 

i<j<k<l 

Thus, we have obtained a class of exact cosmological type solutions for the Gauss- 
Bonnet model (13. ip that is given by the metric (12. 2p with the functions (3 l (u) and 7(w) 
from (13. 8p and (13.91) . respectively, and integration constants c l obeying (I3.10p . 

3.1 Solution with power- law dependence of scale factors 

Let us consider the solutions with 



E c V0. (3.11) 



i=l 



Introducing the synchronous-type variable 



where 



r = - exp(cw + Co), (3-12) 



ri 1 n 

C =gE Ci ' C 0=gE C 0' ( 3 - 13 ) 

i=l i=l 



and defining new parameters 

^ = cVc,, A i = exp[ci+p i (\nc-c )], (3.14) 
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i = 1, . . . , n, we get the "power-law" solution with the metric 

n 

g = wdT®dT + Y^ £ i A i T ^ d y l ® d v\ ( 3 - 15 ) 



1=1 

where w = ±1, £j = ±1; A, > are arbitrary constants, and parameters p J obey the 
relations 

n 

j=i 

Gfy*iPVW = 24 ^ pVpV = 0, (3.17) 

i<j<fc<Z 

following from (I3.10p . (13. lip and (I3.14p . This solution is a singular one for any set of 
parameters p l , see Appendix D. 

In the cosmological case when w = —1, e, = 1 (for all i), this solution was obtained 
earlier in [TU] for D = 5, 6 and verified recently in [2B] for all D > 4. 

Example 1. Let us consider the case D = 6 and pj ^ 0, z = 1, . . . , 5. Relations (13.161) 
and (I3.17P read in this case as follows 

p 1 +p 2 +p 3 +p* +p 5 = 3, (3.18) 

pVpYp 5 ( I + 1 + 1 + 1 + 1 ) = 0. (3.19) 
\p l p z p 6 p° J 

Let us put p 1 = x > 0, p 2 = -, p 3 = z > 0, p 4 = y < 0, p 5 = -. Then we get 

x y 

x + - + z + y + - = 3, x + - + - + y+- = 0. (3.20) 

x y x z y 

Subtracting the second relation in (I3.20p from the first one we obtain 2-| = 3orz = 
|(3 + vT3) (z > 0). For any x > there are two solutions y = y±(x) = ~(— A± ^/ A 2 — 4), 
where A = x + - + ->2. 

X z 

Proposition 1. For D ^ 4 the metric Ii3.15\) is a solution to equations of motion li3.S\) 
if and only if the set of parameters p = (p 1 , ...,p n ) either obeys the relations 113.16}) and 
{3.11^ , or p = (a, b, 0, 0), (a, 0, b, 0, 0), . . where a and b are arbitrary real numbers. 



This proposition is proved in Appendix E. (For cosmological solutions in dimensions 
D = 5,6 see also [ID].) 

For D = 4 the metric (I3.15P gives a solution to equations of motion (13. 2 p for any set 
of parameters p 1 . 

3.2 Solution with exponential dependence of scale factors 

Now we consider the solution with 



8 



J>* = 0. (3.21) 



Introducing the synchronous-type variable 

t = u exp(co), (3.22) 
where cq is defined in (13. 131) and defining new parameters 

v< = d exp(-co), Bi = exp(4), (3.23) 

i — 1, . . . , n, we are led to the cosmo logical- type solution with the metric 

n 

g = wdT®dT + J2 SiB-e 2 "* 1 'dy* <8> dy\ (3.24) 

i=l 

where w = ±1, = ±1; Bi > are arbitrary constants, and parameters v l obey the 
relations 

n 

XV = 0, (3.25) 

1=1 

G» iJH vVW = 24 uVW = 0, (3.26) 

i<j<fc<Z 

following from flBTTUD . (l3T2Tj) and fl3T23|) . 

Example 2. Let .D = 6 and v % ^ 0, i = 1, . . . , 5. Relations ( 13. 25ft and ( 13.261) read in 
this case as follows 

v 1 + v 2 + v 3 + v 4 + v 5 = 0, (3.27) 

„ WvV (1 + 1 + 1 + 1 + 1) =0. (3.28) 
\+ + ir try 

We put v 1 = x > 0, v 2 = -, v 3 = 1, f 4 = y < 0, t> 5 = -. Then we get 

x + - + l + y + - = 0, (3.29) 
For any x > there are two solutions y = y±(x) = \{—B ± y/B 2 — 4), where 5 = 

X + i + 1 > 3. 

x — 

3.3 Some other solutions 

The solutions to equations of motion (13 .6p and (13.71) are not exhausted by relations ( 13. 8ft - 
( 13. 1 0[) . We give an example of another solution for D > 4: 

e-tnr+P+fP+fPpppa = C , (3.30) 
(3\u) = i> 3, (3.31) 
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where /3q (i > 3) and C are arbitrary constants. In terms of "synchronous" variable r 
(obeying dr = e 1<yU ^du) this solutions reads as follows 



g = wdr (g> dr + 2J ^i^1{r)dy l ® dy % 



(3.32) 



where 



dai \ / da 2 \ ( da 3 . 
dr J \ dr J \ dr 







(3.33) 

Oj(Tj = a^, i > 3, (3.34) 

where > (z > 3) and C are constants. This solution contains a special solution with 

<2i(t) = 0.2(1") = as(r) = At. (3.35) 

For C = we get a special solution with arbitrary (smooth) functions 7(w), /3 1 (m), 
fi 2 {u) and constant /^(w) = /3q, for i > 2. In terms of synchronous variable this solution 
is described by the metric (13.321) with 

ai(r), a 2 (r) — arbitrary, aj(r) = a° — constant, i > 2. (3.36) 



Remark 1. For D = 4, or n = 3, the equations of motion (I3.6P and (13. 7p are satisfied 
identically for arbitrary (smooth) functions (3 l (u) and 7(u). This is in agreement with 
that fact that in dimension D — 4, the action (13. ip is a topological invariant and its 
variation is identically zero. 



4 Reduction to an autonomous system of first order 
differential equations 

Now we put 7 = 0, i.e. "the synchronous-like" time gauge is considered. We denote 
u = r. By introducing "Hubble-like" variables h l = we rewrite eqs. (I2.20p and (I2.2ip 
in the following form 



wctiGiMh 1 + a2G ijk ih l h :) h k h l = 



-IwaxGijb? - -a 2 G ijH h J h k h l 



d 

dr 



-■2ic(\ l G ll lv> - ^a 2 G ijkl h j h k h l 



L = 0, 



i = 1, . . . , n, where 



L = -waxGijUb? - -a 2 Gij k ih l h j h k h! . 
3 



(4.1) 



(4.2) 



(4.3) 
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Due to flU]), 

L = —waiGijtfh*. (4.4) 
3 

Thus, we are led to the autonomous system of the first-order differential equations on 
h\r),...,h n {r). 

Here we may use the relations (12.101) . (12. lip and the following formulas (with v % = h l ) 

G i3 ^ = v i -Sx, (4.5) 
G m v j v k v l = Sf + 2S 3 - 33 1 3 2 + 3(S 2 - Sly + 6S 1 {v 1 ) 2 - 6(^) 3 , (4.6) 

i = 1, . . . , n, where S^ = Sk(v) = J2i=i( v *) k - Relation (14. 6p is derived in Appendix B. 

Let us consider the fixed point of the system (14. ip and (14. 2p : h l (r) = v % with constant 
v l corresponding to the solutions 

F = v i T + ? , (4.7) 

where /3q are constants, i = 1, . . . ,n. In this case we obtain the metric (I3.24p with 
exponential dependence of scale factors. For ct\ = we get the solution (I3.24p - (l3.26p . 

Now we put cti 7^ and a 2 ^ 0. For the fixed point v = (v l ) we have the set polynomial 
equations 



G lj v i v j - a w G ijkl v i v j v k v l = 0, (4. 

I 

ci .,(-!, ./.,/"'>"" r' \ r' — ^ 



2GijV J - -a w G ijM v j v k v l 



^-|g^V = 0, (4.9) 

J i=l 



i = 1, ... ,n, where a w = a2{—w)/a\. For n > 3 this is a set of forth-order polynomial 
equations. 

The trivial solution v = (v l ) = (0, 0) corresponds to a flat metric g. 

For any non-trivial solution v we have Y17=t v% 7^ (otherwise one gets from (14.91) 

GyvV = Etli^) 2 ~ (E?=l = aIld heI1Ce V = (°» • • • > ))- 

Let us consider the isotropic case f 1 = ... = t>™ = a. The set of equations (14. 8 j) and 
(14. 9 p is reduced to the equation 

n(n — l)a 2 + a w n{n — l)(n-2)(n-3)a 4 = 0. (4.10) 

For n = 1, a is arbitrary and a = for n = 2, 3. When n > 3, the non-zero solution to 
eq. (I4.10p exists only if a w < and in this case 



^\a w \{n-2){n-3) 



(4.11) 



In cosmological case w = —1, this solution takes place when a^/ai < 0. 

Here the problem of classification of all solutions to eqs. (14. 8p . (14. 9 p for given n arises. 
Some special solutions of the form (a, a, b, b), e.g. in a context of cosmology with 
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two factor spaces, for certain dimensions were considered in literature. See, for example, 

punning. 

Here we outline three properties of the solutions to the set of polynomial equations 

(BP, OH]). 

Proposition 2. For any solution v = (v 1 , ...,v n ) to polynomial eqs. and ^4-9$ : 

i) the vector —v = (—v 1 , — v n ) is also a solution; 

ii) for any permutation a of the set of indices {1, n} the vector v = (v a ^\ ...,v a ^) 
is also a solution; 

iii) there are no more than three different numbers among v 1 , ■■■,v n , whenv = (v 1 , ■■■,v n ) 7^ 
(0,...,0). 

Proof. The first item of the proposition is trivial. The second one follows just from 
relations fl2TTUD . f[2TTTj) . (jPjl and (Q)l . 

Now we prove the item iii). Let us suppose that there exists a non-trivial solution 
v — (v 1 , v n ) with more than three different numbers among v 1 , v n . Due to (14.61) . (14. 9 p 
and Y^i=i v% 7^ an y number v l obeys the cubic equation Cq + C\v l + C2(v l ) 2 + C^{v l )' i = 0, 
with C3 7^ 0, i — 1, . . . , n, and hence at most three numbers among v l may be different. 
Thus, we are led to a contradiction. The proposition is proved. 

This implies that in a future investigations of solutions to eqs. (14 .8p and (14. 9 p for 
arbitrary n we will need a consideration of three non-trivial cases when 1) v = (a, ...,a) 
(see (14. lip ); 2) v — (a, a, b, b) (a 7^ b); and 3) v = (a, a, b, b, c, c) (a 7^ b, 
b 7^ c, a 7^ c). One may put also a > due to item i). 



5 The generalization to the Lowelock model 

The action (II. ip is a special case of the Lowelock model 



S = Jj D z^\^a k C k Y (5.1) 



where a%, a m are constants and C k are defined as follows 

L, k — 4 o Ni N2k n MiM2 ■■■- n M 2k _ 1 M 2k > 

k = 1, . . . , m. (Usually, m is chosen as follows: m = m(D) = [(D — l)/2]; the terms with 
k > m(D) will not give contributions into equations of motion.) Here 

«' = E^' (1) ...^, (5.3) 

cr 

is a generalized Kronecker tensor, totally antisymmetric in both groups of indices: Mi, M 2k 
and Ni, N 2k - In (15.31) a sum on all permutations of the set of indices {1, 2k} is as- 
sumed. Here e a = ±1 is the parity of the permutation a. 

It may be verified that £1 = R[g] and £2 (from (I5.2p ) is coinciding with the Gauss- 
Bonnet term (11.21) . 
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5.1 The Lagrange approach 

Here we suggest the following conjecture: the equations of motion for the Lowelock ac- 
tion (15. ip when the metric (12.21) is substituted are equivalent to the Lagrange equations 
corresponding to the Lagrangian (for cosmological case see also [TOl 128]) 

rn 

L = ^2a k L k , (5.4) 

k=l 

where 

L k = fi k exp[-(2k - 1) 7 + lo]Gt^ A J h • • • ( 5 - 5 ) 
7o = Y^i=i^ l i are rational numbers (/ii = —w,fi2 = —1/3) and 

cl= n (^-!) ( 5 - 6 ) 

l<r<s<2fc 

are the components of Finslerian 2/c-metric: < t>i,...,f 2 fc >2k= Gi*}* k v% i ■■• v% 'tft , v s = 
(v l s ) e R™, where < >2k is a 2fc-linear symmetric form on R™, k = l,...,m. Here 

GfX = G lll2 and = G ili2isU , see fl2ZD and fl23]). 

5.2 Cosmological type solutions for "pure" m-th Lowelock model 

Now we put a>i = ... = a m -\ = and a m ^ 0, i.e. we consider the cosmological type 
model governed by the "pure" m-th Lowelock action 



/ d D zy/\g\Cm,\g] 

J M 



S m = Oi m d zyJ\g\C m [g\, (5.7) 



m — 1, 2, 3, 



It may be verified along a line as it was done in the Section 3 that our conjecture 
implies the existence of cosmological type solutions with the metrics (I3.15P and (I3.24p . 
For the "power-law" solution with the metric 

n 

g = wdr <g) dr + ^ e^T 2 ^ dy l <g> dy i 

i=l 

the parameters p l obey the following relations 

n 

= 2m-l, (5.8) 

i=l 

^l/-^ = ( 2m ) ! Yl P h --P i2m =°- ( 5 - 9 ) 



8l<...<l2n 



instead of (13 . 1 6[) and (I3.17p . (For cosmological solutions see also [TU1 128].) 
For the "exponential" solution with the metric 
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n 

g = wdT®dr + J2 e i Bfe 2v ' T dy i ® dy { 
the parameters v 1 should obey the relations (13.251) : J27=i v% = anc ^ 

Gf™l 2m v ix ...v i2m = (2m)! ^ u* 1 ...-!;**" = 0. (5.10) 

«l<...<i2m 

instead of fl3T26|) . 

The existence of these solutions corresponding to the "pure" Lowelock action (15.71) 
may be considered as test for the validity of the conjecture suggested above. 

6 Conclusions 

Here we have considered the (n + l)-dimensional Einstein- Gauss-Bonnet model. For di- 
agonal cosmological type metrics we have reduced the equation of motion to a set of 
Lagrange equations with the Lagrangian governed by two "minisuperspace" metrics on 
]R n : (i) the pseudo-Euclidean 2-metric (corresponding to the scalar curvature term) and 
(ii) the Finslerian 4- metric (corresponding to the Gauss-Bonnet term). The Finslerian 
4-metric is proportional to n-dimensional Berwald-Moor 4-metric. Thus, we have found 
a rather natural and "legitime" application of n-dimensional Berwald-Moor metric in 
multidimensional gravity with the Gauss-Bonnet term, mth For the case of the "pure" 
Gauss-Bonnet model we have obtained two exact solutions: with power-law and exponen- 
tial dependence of scale factors (w.r.t. "synchronous-like" variable). In the cosmological 
case (with w — — 1, E\ — ... — e n — 1) the first (power- law) solution was obtained earlier 
by N. Deruelle for n = 4, 5 [TU] and verified by A. Toporensky and P. Tretyakov (for 
n — 6, 7) [13] and by S. Pavluchenko (for all n) [28]. See also [37] . 

When the "synchronous-like" time gauge was considered the equations of motion 
were reduced to an autonomous system of first order differential equations. It was 
shown that for any non-trivial solution with the exponential dependence of scale fac- 
tors dj(r) = exp(uV), % = l,...,n, there are no more than three different numbers 
among v 1 , ...,v n (if a± ^ and «2 7^ 0.). This means that the solutions of such type have 
a "restricted" anisotropy. Such solutions may be used for constructing of new cosmo- 
logical solutions, e.g. describing accelerated expansion of our 3-dimensional factor-space 
and small enough variation of the effective gravitational constant. For this approach, see 
[33| [36J and references therein. 

We have also proposed (without a proof) a generalization of the EGB effective 
(cosmological-type) Lagrangian to the Lowelock case (in agreement with [1~0~| [28] for cos- 
mological metrics). According to this conjecture a "pure" Lowelock term of m-th order 
in the action gives a contribution to the effective Lagrangian that contains a Finslerian 
2m-metric. This hypothesis implies the existence of cosmological solutions with power- 
law (see also [101 [28] for cosmological case) and exponential dependence of scale factors 
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for the case of the "pure" Lowelock model of m-th order. A proof of the conjecture men- 
tioned above may be the subject of a separate publication. Another generalization of the 
approach suggested in this paper will be connected with inclusion of a scalar field. 

Here an open problem arises: do the generalized solutions (for arbitrary n) with "jump- 
ing" parameters p l , A4 appear as asymptotical solutions in EGB model when approaching 
a singular point? Recall that Kasner-type solutions with "jumping" parameters p l ,Ai 
describe an approaching to a singular point in certain gravitational models, e.g. with 
matter sources, see [38j EHJ SB El S21 H31 HU H3 H6] and references therein. This prob- 
lem may be a subject of separate investigations. (Here it is worth to mention the paper 
of T. Damour and H. Nicolai [J7], which includes a study of the effect of the 4th order 
in curvature gravity terms, including the Euler-Lovelock term octic in velocities, and its 
compatibility with the Kac- Moody algebra Ei .) 
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Appendix 

A Useful relations for (1 + n)-splitting 

Let us consider the metric defined on K* x MJ 1 (R* = (u-,u+) is an open subset in K) 

n 

g = we 2l{u) du ®du+^ h ij (u)dy i <g> dy j . (A.l) 

Here (hij(uj) is a symmetric non-degenerate matrix for any «6l„ smoothly depen- 
dent upon u. The function 7(u) is smooth. 

The calculations give the following non-zero (identically) components of the Riemann 
tensor 

RoiOj = ~ R%ooj — —Roijo — Riojo — t[ — 2/ijj + 2jhij + hikh kl hij], (A. 2) 

Rijki = ^{-w)e~ 2l (h ik h j i - hiihjk), (A.3) 

i,j,k,l = 1, . . . , n, where here and in what follows h^ 1 = (h^) is the matrix inverse to 
the matrix h = (hij). Here we denote A = dA/du etc. 
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For non-zero (identically) components of the Ricci tensor we get 

#00 = \[-h u hu + l -W'h jk h kl h H + h ik h k a], (A.4) 
Rij = ^(-^)e- 27 [2^- + h i3 {h kl h lk - 27) - 2h lk h kl h l3 ], (A.5) 

i,j = l,...,n. 

The scalar curvature reads 

R = i(-w) e - 27 [4tr(/i/i- 1 ) + tr^/r^tr^/T 1 ) - 4 7 ) - 3tr (hhr^hbr 1 )]. (A.6) 



Let us denote 
(h = (hij)), then 
We obtain 

where 

\h\ = \det(hij)\ and 



M = hh-\ (A.7) 

M + M 2 = hhT 1 . (A.8) 
dfx 



R v ^\ = L 1 + £- } (A.9) 
L x = ^{-w)e-^^/\h\[tiM 2 - (trM) 2 ], (A.10) 



/! = (-w)e-^^/\h\trM. (A. 11) 

In derivation of flA.9j) the following relations were used: 

^M = ^^\h\tT(hh- 1 ), > /H=eVR- (A.12) 
The calculations give us the following relations for quadratic invariants 

RmnpqR MNPQ = -e" 47 {(trM 2 ) 2 - trM 4 + 2tr(2M + M 2 — 2 7 M) 2 }, (A.13) 
8 

RmnR MN = — e" 47 {[-2trM - trM 2 + 2 7 trM] 2 + tr[2M + (trM - 2 7 )M] 2 }. (A.14) 
16 

Relations ( 1A.6[) . (IA. 13[) and (IA.14|) imply the following formula for the Gauss-Bonnet 
term (JOJ 

£2 = — e" 47 {2(trM 2 ) 2 - 2trM 4 + [(trM) 2 - trM 2 ][8trM + 3trM 2 
16 

+ (trM) 2 - 8 7 trM] + 4tr[(M 2 - (trM)M)(4M - 4 7 M + M 2 + (trM)M)]}. (A.15) 
Relation (IA.15I) implies another important formula 

C 2 VW\ = L 2 + -^-f 2 , (A.16) 
16 



where 



and 



L 2 = -Le~ 37 v^{6trM 4 - 3(trM 2 ) 2 
48 

+6trM 2 (trM) 2 - 8(trM)trM 3 - trM 4 } 



/ 2 = ie~ 37 v^{2trM 3 - 3(trM)trM 2 + (trM) 3 }. 



Diagonal metrics. 

Now we consider the diagonal metric 



h^u) = e^eA 



(A. 17) 
(A.18) 

(A.19) 



Si = ±1, i = 1, . . . ,n. Then, = 2(3 l 5ij and we get the following relations for "La- 
grangians" 



-w e 



-7+70 



E(^) 2 -(E^ 



8=1 



i=l 



(A.20) 



1 



i=l i=l j'=l 



n n 



+3(E(^) 2 ) 2 + s(E ^) E(^') 3 - 6 E(^) z 
t=i i=i 



(A.21) 



8=1 



i=l 



where 70 = Yn=iP l - 



The "f-functions" flA.llf) and f|A.18jl read as follows 

n 

A = 2(-u;)e-^^/3 l , 



8=1 



/ 3 = 36 



-37+70 



n n 



2j2(h 3 - KJ2 P) E(^') 2 + E ^ 

8=1 j=i 



i=l 



8=1 



(A.22) 
(A.23) 



B Useful relations for Finslerian 4-metric 



Here we consider a proof of identity ( 12. lip . We decompose the product of 6 terms in the 
definition of the 4-metric ( 12.71) into the sum (of "powers of 5-s" ) 



Gijkl - E G l 



ijkl 



(B.l) 



a=0 



where 



/~i0 1 

^ijkl — L ' ^ijkl 



-5ij — 5ik — dii — 5jk — 6ji — Ski, •••) Gijki ~ ^ij^ikSuSjkSji5ki- 



17 



Then we get 

6 

T = G ijkl vVv k v l = J2 T ^ (B.2) 



a=0 



where T a = G°; jkl v l v j v k v l . 

The calculations of T a give us the following results: 

T° = Sf, T 1 = -6SfS 2 , T 2 = 3Si + 125^3, 

T 3 = -45^3 - 16S 4 , T A = 1554, T 5 = -6S 4 , T 6 = S 4 , (B.3) 

where 

5 fc = s k (v) = ( B - 4 ) 

i=l 

k = 1,2,3,4. 

The summation of all T a in (IB. 3ft leads us to the relation 

T = G iiW vVW = - QS 2 S 2 + 35 2 2 + 8SiS 3 - 6S 4 (B.5) 

coinciding with (12.111) . 

Now we prove relation (I4.6p . We get 

6 

P i = G im v i v k v l = y £ i P?, (B.6) 



a=0 



where P" = G^ jkl v j v k v l , i — 1, . . . 



n. 



The calculations of P? give us the following formulas 

P° = S!, Pl = - 3S 1 S 2 , P 2 = 3£ 3 + 3£ 2 « < + 9S 1 (v i ) 2 , 

P* = -S 3 -3S 1 (v i ) 2 -16(v i ) 3 , if = 15(v < ) 3 , if = -6(V) 3 , ff = (?/) 3 , (B.7) 

z = 1, . . . , n. 

The summation of all P? in (IB.7P leads us to the relation 

Pi = Si + 2S 3 - 3SxS 2 + 3(S 2 - S 2 y + QS^v 1 ) 2 - 6(i/) 3 , (B.8) 

i — 1, . . . , n, coinciding with (14.61) . This relation implies Pfl 1 = T in agreement with the 
definitions (lR2j) and (IB~6l) . 



C Lagrange equations 

Here we prove the relations (12. 1 7[) - (12. 19[) for the cosmological type metric (I2.2p defined 
on manifold M from (12. ip . The tensor £mn is obtained from the variation of the action 



S = [ d D z^/\f\C\g), (c.i; 

J M 
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with C[g] = ai-%] + a 2 C 2 [g], i.e. 



5S= / d D z^\g\S MN 6g MN } (C.2) 

J M 

and v^^a/JV = SS/6g MN . 

Without loss of generality any 1- dimensional submanifold Mj is chosen to be compact 
and coinciding with the circle of unit length: Mj = (r = 1/2tt is the radius of the 
circle) and all coordinates y % (see (12. 1 P ) obey < y l < 1, i — 1, . . . ,n. 

Here we will use the following relations for the components £mn m coordinates (y M ) = 
(y° = u,y l ) and £ calculated for the metric (12. 2ft : 

£-mn = $mn£nn, (C.3) 
£mn = £mn(u), (C.4) 

VW\£ = L + f, (C.5) 
au 

where L = L(j,P,$) and / = f(j,^,$) are defined in relations (12. 4p and (I2.9p . respec- 
tively. 

The first relation (1C.3[) may be readily verified using (I2.14p - (l2.16p and formulas for 
the Riemann tensor (1A.2I) and (1A.3|) . The second relation ( 1C.4P is an obvious one and the 
third one (1C5j) is coinciding with ( 12 .31) . 



The substitution of the metric H 2 . 2 1) into the functional (IC.ip gives us (due to (1C.5P 
and < y l < 1) 

/ df 



S=j du[L + ^-) (C.6) 



and hence 



where 5^f{u) and 5fi l (u) are smooth functions with compact support in (u_, u + ) {5 r y{u±) = 
8(3 l {u±) — 0), i — 1, . . . , n. On the other hand, using flC.2p - flC.4p . the relation 

(^ MiV ) = diag(-2 we - 2 ^ 7 , -2e 1 e- 2/3l (5/3 1 , . . . , -le^ 6f3 n ) 
and < y l < 1, we get 

rfu{^ 00 (-2w)e 7 °- 7 57 + ^ ^(-2e i )e 7+7 °- 2/3I 5/3 4 }. (C.8) 

i=i 

Comparing flUTj) and (ICT8]l we get relations f l2TTj) and fl2TT8|) . Relations fl2TT9|) just 
follow from f lU3|) . 
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D Riemann tensor squared 

Here we consider the Riemann tensor squared (Kretchmnann scalar) for the metric (I3.15P 

n 

g = wdr ®dr + ^2 s i A 2 i T 2pl dy l <g> dy\ 
i=i 

From (1A.13I) we get 

Rmnp Q R MNPQ = Kr~\ (D.l) 

where 

= 25 4 + 25f - 85 3 + 45 2 (D.2) 

and S k = S k (p) = ElUW* k = 1,2,3,4. 
Using the identities 

n 

K = 4 - lf(ff + 2(5 2 2 - 5 4 ). (D.3) 

i=i 

and 

5 2 2 - 5 4 = 2^>*)V) 2 (D.4) 

we obtain that X > and = if and only if the set of parameters p = (p l , ...,p n ) is 
either trivial: p — (0, 0), or belongs to the Milne set: 

p=(l,0,...,0),...,(0,...,0,l). (D.5) 

For other sets p we have K > and the Riemann tensor squared diverges when 
r +0. 



E The proof of Proposition 1 

The equations of motion (14. ip and (14. 2 j) corresponding to the metric (I3.15P with /i J = p l /r 
(here «i = and «2 7^ 0) read as follows 

A = G ijklP yp k p l = Q, (E.l) 



A = G ljkl p>p k p l = 0, (E.2) 



z = 1, . . . , n. 

Let D = n + l/ 4 and 



n 



1 

3< 



£=1(8-2}), (E.4) 
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i = 1, . . . ,n. 

For D 7^ 4 the set of equations (lE.ip and (1E.2[) is equivalent to the following set of 
equations 

A = St - 6S 2 S 2 + 3S 2 2 + 8S1S3 - 6S 4 = 24 ^ pVpV = 0, (E.5) 

i<j<fc<Z 

B = (5i - 3)(Sf - 3SiS 2 + 2S 3 ) = 6(S 1 - 3) ^ pVp* = 0, (E.6) 

«<i<fc 

& = (S 1 - 3)pW) 2 - 2S lP i + Sf - S 2 ) = 0, (E.7) 

i = 1, . . . ,n. Here 5*^ = Skip) = Y^=i(P z ) k aric ^ we use d the identities (12.111) . (14. 6 p and 
the following identity 

Sf - 3S.S2 + 2S 3 = G ijkP y P k = 6 (E.8) 

i<j<k 

where 

G^ = (^-l)(^-l)(^-l) (E.9) 

are components of a Finslerian 3-metric. The identity flE.8j) could be readily verified along 
a line as it was done in Appendix B for the Finslerian 4-metric. (We note that relation 
( IE. 61) may be also obtained using the formula ( 1A.16I) .) 

For S\ = 3 we obtain the main solution governed by relations (I3.16P and (I3.17p . 

Now we consider another case Si ^ 3. Let k be the number of all nonzero numbers 
among p 1 ,...,p n . For k = we get a trivial solution (0, ...,0). Let k > 1. We suppose 
without loss of generality that p 1 ,...,p k are nonzero. For k = 1,2 all relations (IE. 5ft - 
(IE.7I) are satisfied identically. In all three cases k — 0,1,2 the solutions have the form 
(a, b, 0.., 0) (plus permutations for general setup). 

Now we consider k > 3. From ( IE.7j) and Si ^ 3 we obtain 



2(p l Y - 2S lP l + Sf - S 2 = 0, (E.10) 

i — 1, . . . , k. Summing on % gives us (2 — k){S 2 — Sf) = 0, or S2 = Sf. Then we obtain 
from (IE.6P S3 = Sf and from f)E.5[) : S4 = Sf. Thus, we get S 4 = S| implying S = 
Xa^^xfcO 9 *) 2 ^) 2 = 0- S > (p x ) 2 (p 2 ) 2 > 0. Hence, we are led to a contradiction. 
That means that for Si 7^ 3, we have only solutions with k < 2 of the form (a, b, 0..,0) 
(plus permutations for general setup). The Proposition 1 is proved. 

References 

[1] B. Zwiebach, Curvature squared terms and string theories, Phys. Lett. B 156, 315 (1985). 

[2] D. Gross and E. Witten, Superstrings modifications of Einsteins equations, Nucl. Phys. B 
277, 1 (1986). 



21 



[3] D.J. Gross and J. H. Sloan, The quartic effective action for the heterotic string, Nucl. Phys. 
B 291, 41 (1987). 

[4] R.R. Metsaev and A. A. Tseytlin, Two loop beta function for the generalized bosonic sigma 
model, Phys. Lett. B 191, 354 (1987). 

[5] R.R. Metsaev and A. A. Tseytlin, Order alpha-prime (two loop) equivalence of the string 
equations of motion and the sigma model Weyl invariance conditions: dependence on the 
dilaton and the antisymmetric tensor, Nucl. Phys. B 293, 385 (1987). 

[6] H.-J. Schmidt, Fourth order gravity: Equations, history, and applications to cosmology, 
Int. J. Geom. Meth. Mod. Phys. 4, 209-248 (2007); arXiv: gr-qc/0602017] 

[7] S. Nojiri and S.D. Odintsov, Introduction to modified gravity and gravitational alternative 
for Dark Energy, Int. J. Geom. Meth. Mod. Phys. 4, 115-146 (2007); |hep-th/0601213| 

[8] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov and S. Zerbini, One-loop effective action 
for non-local modified Gauss-Bonnet gravity in de Sitter space, arXiv: 0905.0543. 

[9] H. Ishihara, Cosmological solutions of the extended Einstein gravity with the Gauss-Bonnet 
term, Phys. Lett. B 179, 217 (1986). 

[10] N. Deruelle, On the approach to the cosmological singularity in quadratic theories of gravity: 
the Kasner regimes, Nucl. Phys. B 327, 253-266 (1989). 

[11] E. Elizalde, A.N. Makarenko , V.V. Obukhov , K.E. Osetrin and A.E. Filippov, Stationary 
vs. singular points in an accelerating FRW cosmology derived from six-dimensional Einstein- 
Gauss-Bonnet gravity, Phys. Lett. B 644, 1-6 (2007); |hep-th/0611213| 

[12] K. Bamba, Z.-K. Guo and N. Ohta, Accelerating Cosmologies in the Einstein-Gauss-Bonnet 
theory with dilaton, Prog. Theor. Phys. 118, 879-892 (2007); arXiv: 0707.4334. 

[13] A. Toporensky and P. Tretyakov, Power-law anisotropic cosmological solution in 5+1 di- 
mensional Gauss-Bonnet gravity, Grav. Cosmol. 13, 207-210 (2007); arXiv: 0705.1346. 

[14] I.V. Kirnos, A.N. Makarenko, S.A. Pavluchenko and A.V. Toporensky, The nature of singu- 
larity in multidimensional anisotropic Gauss-Bonnet cosmology with a perfect fluid, arXiv: 
0906.0140. 

[15] S.A. Pavluchenko and A.V. Toporensky, A note on differences between (4+1)- and (5+1)- 
dimensional anisotropic cosmology in the presence of the Gauss-Bonnet term, Mod. Phys. 
Lett. A 24, 513-521 (2009). 

[16] I.V. Kirnos and A.N. Makarenko, Accelerating cosmologies in Lovelock gravity with dilaton, 
arXiv: 0903.0083. 

[17] M. Kowalski, D. Rubin et al., Improved cosmological constraints from new, old and com- 
bined supernova datasets, arXiv: 0804.4142. 

[18] D.G. Boulware and S. Deser, String generated gravity models, Phys. Rev. Lett. 55, 2656 
(1985). 



22 



J.T. Wheeler, Symmetric solutions to the Gauss-Bonnet extended Einstein equations, Nucl. 
Phys. B 268, 737 (1986). 

J.T. Wheeler, Symmetric solutions to the maximally Gauss-Bonnet extended Einstein equa- 
tions, Nucl. Phys. B 273, 732 (1986). 

D.L. Wiltshire, Spherically symmetric solutions of Einstein-Maxwell theory with a Gauss- 
Bonnet term, Phys. Lett. B 169, 36 (1986). 

R.-G. Cai, Gauss-Bonnet black holes in AdS spaces, Phys. Rev. D 65, 084014 (2002). 

M. Cvetic, S. Nojiri and S.D. Odintsov, Black hole thermodynamics and negative entropy in 
de Sitter and anti-de Sitter Einstein-Gauss-Bonnet gravity, Nucl. Phys. B 628, 295 (2002). 

C. Garraffo and G. Giribet, The Lovelock black holes, Mod. Phys. Lett. A 23, 1801 (2008); 
arXiv: 0805.3575. 

C. Charmousis, Higher order gravity theories and their black hole solutions, Lect. Notes 
Phys. 769, 299 (2009); arXiv: 0805.0568. 

K.A. Bronnikov, S.A. Kononogov and V.N. Melnikov, Brane world corrections to Newton's 



law, Gen. Rel. Grav. 38, 1215-1232 (2006); arXiv: [gr^qc/ 0601114. 



E. Kasner, Geometrical theorems on Einstein's cosmological equations, Amer. J. Math., 
43, 217 (1921). 

S.A. Pavluchenko, On the general features of Bianchi-I cosmological models in Lovelock 
gravity, Phys. Rev. D 80 (2009), 107501; arXiv: 0906.0141. 

L. Berwald, Projective Krummung allgemeiner affiner Raume und Finslersche Raume 
skalarer Krummung, Ann. Math. 48, 755-781 (1947). 

A. Moor, Erganzung, Acta Math. 91, 187-188 (1954). 

G. Bogoslovsky, Rapidities and observable 3- velocities in the flat finslerian event space with 
entirely broken 3D isotropy, SIGMA 4, 045, 21 pp. (2008); arXiv: 0712.1718. 

G.I. Garas'ko and D.G. Pavlov, Construction of the pseudo Riemannian geometry on the 
base of the Berwald-Moor geometry; |math-ph /0609009 , 

V.D. Ivashchuk and V.N. Melnikov, Exact solutions in multidimensional gravity with 
antisymmetric forms, topical review, Class. Quantum Grav. 18, R82-R157 (2001); 
|hep-th/0110274| 

D. Lowelock, The Einstein tensor and its generalizations, J. Math. Phys., 12, No. 3, 498-501 
(1971). 

V. Baukh and A. Zhuk, Sp-brane accelerating cosmologies, Phys Rev. D 73, 104016, 21pp. 
(2006). 

V.D. Ivashchuk, S.A. Kononogov and V.N. Melnikov, Electric S-brane solutions correspond- 
ing to rank-2 Lie algebras: acceleration and small variation of G , Grav. Cosmol. 14, No. 
3, 235-240 (2008); arXiv: 0901.0025. 



23 



[37] V.D. Ivashchuk, On anisotropic Gauss-Bonnet cosmologies in (n + 1) dimensions, governed 
by an n-dimensional Finslerian 4-metric, arXiv: 0909.5462. 



[38] V.A. Belinskii, E.M. Lifshitz and I.M. Khalatnikov, Oscillating regime of approaching to 
peculiar point in relyativistic cosmology, Usp. Fix. Nauk 102, 463-500 (1970) (in Russian); 
A general solution of the Einstein equations with a time singularity, Adv. Phys. 31, 639-667 
(1982). 

[39] J. Demaret, M. Henneaux and P. Spindel, No oscillatory behavior in vacuum KaluzaKlein 
cosmologies, Phys. Lett. 164 B, 27-30 (1985). 

[40] V.D. Ivashchuk, A. A. Kirillov and V.N. Melnikov, On stochastic properties of multidimen- 
sional cosmological models near the singular point, Izv. Vuzov (Fizika) 11, 107-111 (1994) 
(in Russian) [Russian Physics Journal 37, 1102 (1994)]. 

[41] V.D. Ivashchuk, A. A. Kirillov and V.N. Melnikov, On stochastic behaviour of multidimen- 
sional cosmological models near the singularity, Pis'ma ZhETF 60, No 4, 225-229 (1994) 
(in Russian) [JETP Lett. 60, 235 (1994)]. 

[42] V.D. Ivashchuk and V.N. Melnikov, Billiard representation for multidimensional cosmology 
with multicomponent perfect fluid near the singularity, Class. Quantum Grav. 12, No 3, 
809-826 (1995); |gr^qc/9407028[ 

[43] V.D. Ivashchuk and V.N. Melnikov, Billiard representation for multidimensional cosmology 
with intersecting p-branes near the singularity, J. Math. Phys. 41, No 9, 6341-6363 (2000); 
|hep-th/9904077| 

[44] T. Damour and M. Henneaux, Chaos in superstring cosmology, Phys. Rev. Lett. 85, 920-923 
(2000); |hep-th/0003139| 

[45] T. Damour, M. Henneaux and H. Nicolai, Cosmological billiards, topical review, Class. 
Quantum Grav. 20, R145-R200 (2003); [hep^th/0212256[ 

[46] V.D. Ivashchuk and V.N. Melnikov, On billiard approach in multidimensional cosmological 
models, Grav. Cosmol. 15, No. 1, 49-58 (2009); arXiv: 0811.2786. 

[47] T. Damour and H. Nicolai, Higher order M theory corrections and the Kac-Moody algebra 
E10, Class. Quantum Grav. 22, 2849-2879 (2005). 



24 



